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We investigate the properties of forced inertial modes of a rotating fluid inside a spherical 
shell. Our forcing is tidal like, but its main property is that it is on the large scales. By 
numerically solving the linear equations of this problem, including v iscosity, we fir st con- 
firm some analytical results obtained on a two-dimensional model bv lOgilviel(|200,^ : some 



additional properties of this model are uncovered like the existence of narrow resonances 
associated with periodic orbits of characteristics. We also note that as the frequency of 
the forcing varies, the dissipation varies drastically if the Ekman number E is low (as 
is usually the case). We then investigate the three-dimensional case and compare the 
results to the foregoing model. The three-dimensional solutions show, like their 2D coun- 
terpart, a spiky dissipation curve when the frequency of the forcing is varied; they also 
display small frequency intervals where the viscous dissipation is independent of viscosity. 
However, we show that the response of the fluid in these frequency intervals is crucially 
dominated by the shear layer that is emitted at the critical latitude on the inner sphere. 
The asymptotic regime, where the dissipation is independent of the viscosity, is reached 
when an attractor has been excited by this shear layer. This property is not shared by 
the two-dimensional model where shear layers around attractors are independent of those 
emitted at the critical latitude. Finally, resonances of the three-dimensional model corre- 
spond to some selected least-damped eigenmodes. Unlike their two-dimensional counter 
parts these modes are not associated with simple attractors; instead, they show up in 
frequency intervals with weakly contracting webs of characteristics. Besides, we show 
that the inner core is negligible when its relative radius is less than the critical value 
OAE^^^. For these spherical shells, the full sphere solutions give a good approximation 
of the flows. 



1. Introduction 

1.1. The background 

The question of how close binary stars born on eccentric orbits reach a circular trajectory 
with synchronous orbital and spin rotation is by far not fully resolved. The process of 
circularization of the orbits and synchronization of rotations is thought to be the result 
of the tidal interaction. Tides act in two ways: they first create a tidal bulge which 
generates a torque leading to angular momentum exchange between the orbital motion 
and the spinning stars. Second, they dissipate energy through the induced fluid flows. 
This latter effect is important as it controls the time scale of the whole evolution process. 
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Thus, much work has been devoted to evaluate the efficiency of the dissipative pro- 
cesses, which can be triggered by tides. The first process that has been investigated is 
the viscous damping of the equilibrium tide; this is the tide induced by the mere defor- 
mation of equipotentials; the fluid is assumed "attached" to these su rfaces and the tim e 
dependence of the tidal distortion makes the flow (see ,Zahn 19661 Il977l . 1992 . 2008 ). 
However, this process is efficient if the viscosity is high enough, a condition which is met 
only by low-mass stars. Indeed, the convective envelope of these stars provides a strong 
damping through their turbulent viscosity. However, stars of mass larger than 1.8 solar 
mass have no convection in their outer layers (only in the central part). Their envelope is 
stably stratified. But observations do show that circularization and synchronization are 
also effective for these stars when they are close enough (Giuricin et a/,, 1984 . but see also 
the recent review of observationa l facts on tidal interactions between stars or between 
stars and planets bv iMazetil 20081 ). The present understanding of these results is based 
on the idea that the radiative envelope of these stars are efficiently d amping the g ravity 
waves excited by the tides. This mechanism was first investigated bv IZahn ( 1975f ) using 
a simplified approach with an asymptotic description of gravity modes in a non-rotating 
fluid. However, stars in a binary system are rotating, some times quite rapidly. Rotation 
thus appeared as an unavoidable feature. Following work of lRoccal(ll987lll989l) included 
therefore the Co riolis acceleration but as a first order perturbation. 

More recentlv. I Witte fc Savonii^ (1999a, 2001) investigated numerically the synchro- 
nization process of massive stars, fully taking into account the Coriolis force in the flow. 
They thus discovered the phenomenon of "resonance locking" by which the tidal forc- 
ing excites resonantly two rotationally modified gravity modes (hereafter called gravito- 
inertial modes), which therefore strongly dissipate energy; the locking occurs because one 
mode tends to spin the star up while the other tends to spin it down. The low frequen- 
cies of the resonant gravito-inertial modes, however, is synonymous of short-wavelength 
modes that are not well s uited to the forego i ng two -dimensional numerical calculations. 
This difficulty motivated ISavonije fc Wittel (|2002h to further investigate the problem 
with the so-called Traditi onal Approxima t ion w hich allows a separation of the variables. 
However, as discussed bv lGerkema et al. (2008), this simplificati on eliminates a crucia l 
feature of these modes, namely their internal shear layers. Indeed, Dintrans et al. ( 19991 ) 
showed that gravito-inertial modes are singular in the limit of vanishing diffusivities. 
This property is actually shared by gravity mo des and inertial mod es in any configu- 
ration where variables cannot be separated (see iRieutord et al. 2000[ ) . The singularities 
come from the ill-posed nature of the mathematical problem which is of hyperbolic or 
mixed type with boundary conditions. At non-zero but low diffusivities, dominant singu- 
larities show up as modes confined in shear layers f ollowing attractors of characteristics 
(jRieutord fc V aldettaro 1997; Rieu tord et aLll200lh . As far as dissipation is concerned, 
such modes behave very differently compared to the regular ones. Astrophysical situ- 
ations being characterised by large ratios between integral and dissipation scales, it is 
crucial to understand the asymptotic properties of resonant-mode dissipation at vanish- 
ing viscosities and thermal diffusion. 

Besides the problem of tidal interactions between stars, which is quite old, the recent 
discovery of many planetary systems harbouring Jupiter-like planets on orbits very close 
to the central star also motivates a new examination of the tidal interaction. The nov- 
elty with planets, even of the size of Jupiter, is that they likely contain a rocky core in 
their central part. Even if thi s core is not a solid body, th e transition with the surround- 
ing envelope is likely sharp ( Goodman fc Lackner 20091 ). This makes the fluid domain 
like a spherical shell. In stars this domain is delineated by the Brunt - Vaisala frequency 
variations and may be more complex (e.g. Dintrans fc RieutordI 2000[ l. 
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Beyond the astrophysical problem described above, the case of resonant i nertial modes 



is als o of interest in Earth sciences in rela tion to the elliptic instability (e.g. lLacaze et 



200i), the rotating precessing flows (e.g. iHoUerbach fc KerswelllllQQSf) or for the under 

standing of the dynamics of the ocean or of the atmosphere (e.g. Maasl2001 : Maas &: Harlander 



[2007 '). Furthermore, it has also been investigated in the context of en gineering a,pplica - 
tions such as the oscillations of fuel tanks of spinning spacecrafts (see [Manassehlll996 ). 
Finally, let us mention that non-axisymmetric inertial modes also appear as the growinj 
pertu rbations of a rotating fluid destabilized by thermal convection (see ;Zhana il99 



199i). 



As may be guessed, the full astrophysical problem is very involved and some simplifi- 
cations are in order if one wishes to decipher the mechanisms controll ing the asymp totic 
laws of dissipation at small diffusivities. Thus, following the work of Ogilvie (2003), we 
first investigate the case of forced singular modes in a sle nder toroidal sh ell, which is 
a two-dimensional approximation for a spheric al shell ( e .g. Rieutord et aLi , 20021. Thus 
doing, we can recover and extend the results of [Ogilvid ( 2005 ). We proceed then by fo- 
cusing on the thr ee-dimensi onal problem o f the spherical shell and extend the previous 
work of iRieutordI (,1991) and|Tilgner| (|l999[ ). We then discuss and compare the results of 
2D and 3D models. Some conclusions end the paper. 

1.2. The model 

We consider a viscous fluid inside a spherical shell that mimics a stellar or planetary 
envelope, submitted to the tides of an orbiting mass point. This tidal forcing may be 
condensed in the tidal potential which we write as: 



^P2(C0S6') COS UJot + ^nxf^ff (C0S6') COs{2LUot - 2lf>) 



following IZahnI (|1977[ ). In this expression, uJq is the orbital angular velocity of the point 
mass, (r, 0, f) are the spherical coordinates whose origin is at the centre of the body un- 
der consideration; P2 and P| are Legendre polynomials, while <I>ax and $nx are the ampli- 
tudes. On Earth <I>nx is the amplitude of the well-known semi-diurnal tide. 
$ax'''^P2(cos6') cos Ldot \s the leading term coming from the eccentricity of the orbit. Al- 
though <I>ax is usually smaller than $nx, we shall discard $nx and only consider the first 
term of $t- This is justified by the fact that the properties of non-axisymmetric iner- 
tial modes in a spherical shel l are the same as tho se of the axisymmetric ones as far 
as singularities are concerned ([Rieutord et aLll200ll) . Thus all the results derived below 
for an axisymmetric forcing can be applied, mutatis mutandis, to a non-axisymmetric 
forcing. In another probl em like the res onant interaction of inertial modes, which leads 
to the elliptic instability (|Kerswel]||2002l ) , the non-axisymmetric terms would be essential 
of course. 

The fluid inside the spherical shell is rotating at the spin angular velocity of the 
star. Thus, in a frame co-rotating with the fluid, the tidal forcing is 

$T = «'axr-^P2(c0s6l) COs[(Wo - n)t\ 

The first response of a star to the tidal potential is the so-called equilibrium tide, which 
describes the distortion of the equipotentials. This induces a velocity field Ve, which is 
derived from the time-evolution of the equipotentials. Following ZahnI ( 1966( ). this induces 
a radial flow of the form: 

Ve = Ar^ P2{cos9) sinutcr 

where oj = ujo — As shown by Ogilvie ( 20051 ). this flow forces a dynamical response of 
the star through the body force 
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f = -{iLjVe + 2n X Ve) 



We are interested in the dynamical response of the fluid. As the general realistic case 
is much involved, we reduce it to the study of the response of an incompressible viscous 
fluid inside a spherical shell. Although much simplified, this model retains the essential 
feature of the low-frequency ste llar modes, namely their s ingularities associated with 
attractors of characteristics (e.g. Dintrans &: Rieutord 2000l ). 

Assuming that the fluid response remains of small amplitude, we need to solve the 
linearised equations governing forced periodic perturbations of a viscous rotating fluid 
with constant density. When the length scale is the outer radius of the shell R, the 
time scale is (2fl)^^, the equations of the non-dimensional pressure {p) and velocity 
perturbation (m) may be written: 



iuju + Cz y< u 



V • u = 



-Vp + EAu + f 



(1.1) 



where E = v/2VLR^ is the Ekman number and v is the kinematic viscosity. The non- 
dimensional force reads 



f ^ ~Veir,0){iu}er + smee^) (1.2) 

where Ve{r^O) — Ar^P2{cos9). We complete these equations with stress-free boundary 
conditions. The use of stress-free boundary conditions may be surprising in the plan- 
etary case as the interface between the fluid layer and the core is a solid boundary. 
However, because of the nature of the fluid flows that are restricted to internal shear 
layers, this boundary condit i on ha s little effect on the solution as was actually found by 
Fotheringham fc HoUerbachl (|l998r ). We give in appendix B the scaling arguments that 
lead to this result. 

Even thus simplifie d, our problem remains challenging. A further simplific ation, actu- 
ally used many times ( Maas fc Lam 1995 ; Rieutord et al. 2002 ; Ogilvidl2005f ) , consists in 
reducing the problem to two dimensio ns. In such a case, th e spherical shell turns into a 
(cored) slender torus (see figure [T] and iRieutord et a/.ll2002l) . but the eige nvalue problem 



as well as the forced one (far from resonances) can be solved analytically ([Rieutord et 



l2002t lOgilvie 20051 ). Although this latter simplification seems quite drastic, it may be 
shown that it is nevertheless relevant to equatorial regions of a thin sph erical sh e ll (e.g . 
Stewart sonlll971 : Rieutord et "aLl l2Q02l) . To bridge the gap with the work of| iQgilvid (l2005h . 
we therefore first consider this two-dimensional model and then move on to the three- 
dimensional one. 



2. Resonances of axisymmetric modes in the slender torus 

2.1. Mathematical formulation 

Using cylindrical coordinates (s, z), we first introduce a meridional stream function -0 
so that mass conservation is automatically insured. Hence, we set 

dip 1 dsip 

oz s OS 

The momentum equation leads to 
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Figure 1. A schematic view of the slender cored torus: the principal radius a of the torus 

tends to infinity. 



where 



and 



EA'A'i} - d,u - iujA'i) = C 
EA'u + dzTp — ioju = — / 



ds \s ds I dz'^ 



(2.1) 



fs = -iLo sin eVe (r, e) h = -iLU COS OVe (r, 9) 
The torque density C may be further reduced into 

_ iuj dVe 

~ r de 

System p.ip is completed by stress-free boundary conditions which demand that 



d /ldij\ du ^ 
or \r or J or 

Here 77 is the non-dimensional radius of the inner core. 

2.2. Solving the equations 

Because of the shape of the boundaries, it is natural to use polar coordinates (p, 0) in 
a meridional section of the torus (see the schematic view in figure [1]) . These coordinates 
are related to the foregoing cylindrical coordinates by 

s — pcoscf), z = psincj) 

They are also related to the original spherical coordinates by r = p and 9 = ^ — 4>. 

The two-dimensional approximation is essentially summarised into the neglect of the 
curvature terms so that, for instance, tp/s <C dtp/ds; hence. 

Up — Us COS (p + Uz sm ffl = — - 

p ocj) 



COS ( 
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dip 

'dp 

To obtain a numerical solution of the equations (j2.ip (simplified with the 2D-approximation) , 
we first use a Fourier decomposition, namely 

(V-, U, f, C){p, 0) = ^(Vn, -iVr,, fn, C„)(p)e^"*, 



and find the set of ordinary differential equations controlling the shape of the radial 
functions ipn{p) and Vn{p). They are 



2 
V' 



where 



2p 

(Ti-i)y„_i+(«+i)y„+i 

2p 



^ 9^ 1 9 71^ 
" p dp p^ 



With this formulation, the stress-free boundary conditions read: 

d^lpn 1 dlpn _ dVn 

dp 



dp^ p dp 
For the tidal forcing at hands, we find that 







(2.2) 



so that 



C(p,0) 



fn — /— r 



(cos (/) — 2 sin^ (f) cos 0) = ^ (cos </> + 3 cos 3(/)) 
: — 3«a;/9sin0cos0 



— (<5„,i + 3(5„,3), 
16 



C'n 



(2.3) 



The two-dimensional set-u p restricts the iiumb er of excited attractors compared to the 
three-dimensional one fsee iRieutord et a/.ll2002i and below). Thus, although this forcing 
is derived from the true tidal force, we complete our view of the solutions by the use of 
the following forcing 



thus 



/(p,(/.) = 2cos2</), C = 0; 



fn — f-n — Sn^2, Cri — 



2.3. Symmetries 

In modelling the tidal interaction, it is commonly assumed that spin and orbital angular 
momentum vectors are parallel. With this assumption, the tidal flow is symmetric with 
respect to the equatorial and orbital planes. In this case, the tidal force is such that 
/„ — f-n and C„ = —C-n as noted above. 

This symmetry allows us to solve the set of equations (j2.2p solely for n > 0. However, 
solutions in the torus verify a further symmetry: they may be symmetric or antisymmetric 
with respect to the transformation (f) (j) A- it . Such a symmet r y is sp ecific to the torus 
and does not exist in the sphere. As noticed bv lRieutord et al. ( 2002 ). it has a selection 
effect on attractors: some may exist in the spherical shell but be not authorised in the 



Viscous dissipation by tidally forced inertial modes in a rotating spherical shell 7 



torus. For this reason, some attractors of the spherical sheU cannot be studied in two 
dimensions with the natural forcing p.3|) . but can be investigated with the second forcing. 
Still some others cannot be studied in two dimensions altogether. 

The equatorial symmetry of the tidal forcing in combination with the parity of its 
Fourier component, implies that only half of the Fourier components are excited, namely 

T/l, T^3, . . . , V2n+1, ■■■ "02, • • ■ , 1p2n ■ ■ ■ 

we also note that, since V-i = Vi, 

AAtpo ~ = 

which means that ipQ is not excited and therefore vanishes. The second forcing excites 
the other set of Fourier components, namely 

Vo,V2, ■ ■ ■,V2n, ■■■ '01, • • ■ ,'02n+l • ■ ■ 

2.4. Dissipation 

As discussed in the introduction, total viscous dissipation of the fluid volume is the 
actual quantity to be evaluated when the solution is kn own (th i s is th e quantity which 
controls the secular evolution of the orbit). As noted bv lOgilvie ( 2005 ). dissipation may 
be evaluated in two ways: 



D = ^ f = / Re{u* ■ f)dV (2.4) 

2 ./(V) 



where [c] is the rate-of-strain tensor. We use both of these expressions to evaluate the 
internal numerical precision of our results. 

Let us mention that the components of the rate-of-strain tensor [c] in the (p, 4>) coordi- 
nates are easily obtained if these coordinates are completed by an axial one, let say (, so 
that (p, 0, C) form a cylindrical system of coordinates (in our problem, C is perpendicular 
to the meridional plane) ; we show in appendix the equivalence of these expressions with 
the usual ones. In these coordinates, the rate-of-strain tensor components are 



2 - 



1 d^ij 1 dtp\ d^ij) I dip 1 d^ip 



du 1 du 

and the volumic dissipation reads 

D = 2{\cpX + Kc? + M + M) 
2.5. Results 

For comparison with previous work, in all the numerical applications, we set the ratio of 
the inner radius to the outer one r] to 0.35 which is the value of the Earth's liquid core. 

2.5.1. Overall properties 

As shown bv IWitte fc Savoniie (19996), during the evolution of the orbits, the tidal 



forcing frequency scans the whole inertial band and resonances there play a major part. 
Let us point out that in the text below, we shall understand the term "resonance" as the 
local maximum (in frequency) of the viscous dissipation. This definition is appropriate 
for our problem but slightly different from the usual one which refers to the amplitude 
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Figure 2. Scan of the dissipation as a function of the critical latitude angle at _E = 10 ^. 

Resolution is Nr=180, N0=4OO. 
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Second Forcing, Fkman=rO" 
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Figure 3. Same as in figure [2j but for the second forcing. 



of the flow. Resonances are usually understood as the signature of the excitation of 
eigenmodes. As will be clear later, this classical view is not always appropriate here 
because of the ill-posed nature of the inviscid problem. Thus our wording "resonance" 
should be understood in a rather loose way. 

With the tidal forcing in mind, we simulate a scan of the inertial band. We represent 
the viscous dissipation as a function of the frequency, or, more appropriately, as a function 
of the critical latitude i? = arcsin(ci;). As shown by figures [2] and [3l the curve D{uj) is 
very spiky and almost symmetric with respect to 7r/4. This underlying symmetry is a 
consequence of an invariance of Poincare equation in two-dimensions. Indeed, setting 

= in (|2.ip and rewriting the equation for -0, we find 

in the two-dimensional limit. If _F — we easily see that this equation is invariant with 
respect to the transformation s — s- z, z — > s and oj \/l — uP- . Thus the place of resonant 
frequencies is indeed symmetric with respect to 7r/4, if the critical latitude is used as a 
variable. However, the final curves, which are shown in figures [2] and [3l partially lose 
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Figure 4. The attractor occupying the frequency band [0.60922,0.62276]. Ahhough prominent 
in the 3D problem, this attractor cannot be excited by the natural tidal forcing; fortunately it 
can be studied with the second forcing. The small segment indicates the place where the velocity 
profiles of figure [6] have been taken. On the right we show the kinetic energy distribution in a 
meridional plane for a numerical solution computed aX uj — 0.621 and E = 10~^. 




Figure 5. Left; Dissipation as a function of forcing frequency uo near the attractor of figure 
[l] Note that far from the resonance the dissipation tends to be independent of the Ekman 
number E. Right: For the same attractor, scaled resonance curves. Note that when scaled as 
indicated, dissipation curves no longer depend on viscosity. Here uo — ujt is the upper bound of 
the frequency interval of the attractor. 



their symmetry because the forcing does not verify the foregoing invariance as well as 
the viscous terms (the second forcing breaks more strongly this symmetry). 

2.5.2. Attractors and associated resonances 

To better understand the properties of the dissipation curves as shown in figures [2] 
and[3l we shall first concentrate on a given attractor, which is displayed in figure ID This 

attractor may be found in the frequency range [wajt^^b] where ujh = y^ ^-v^"^^ ' 

no analytic expression). In figure [5] we give a detailed view of the resonance associated 

with this attractor around its asymptotic frequency Wf,. 

The remarkable property of these curves is the independence of dissipation from viscos- 
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E = IO-' 
E = IO-" 
E = 10-» 




Figure 6. Profiles of the derivative along the straight thick line shown in figure lUleft 
crossing a segment of the attractor aX u = 0.62f , rj — 0.35. The abscissa is normalised by E^^^: 
note that the location of the maxima of is nearly the same for the three profiles, showing 
that the width of the layer scales with E^^^ . 



ity for frequencies far from th e asym ptotic frequency uji . This result perfectly illustrates 
the demonstration of lOgilvid (|2005l ) who showed that dissipation by inertial modes as- 
sociated with attractors with a finite Lyapunov exponent is actually independent of the 
Ekman number. This comes from the width of the associated shear layers, which scale 
like E^/^. Fi gure [6] indeed confirms that velocity gradients inside a shear layer vary on a 
scale changing with the on e-third power of the viscosity. 

In lRieutord aZ.I (|2002} ) however, it was shown that shear layers associated with freely 
decaying inertial modes, shaped by attractors, have a width scaling like E^/'^. This means 
that at resonances, dissipation resulting from the excitation of such singular modes di- 
verges as E~'^/'^ in the asymptotic limit of sm all Ekman numbers (see below sect. 12.5.^ . 
Moreover, as shown bv iRieutord et al. ( 2002 ). the frequency of these eigenmodes is of 
the form uj = lvq+i^iE^/'^ for <C 1. We thus expect resonance widths to be of the order 
of E^^^. Figure [SJd perfectly illustrates these properties, showing that if scaled properly, 
resonance curves are independent of viscosity. 

2.5.3. Resonances associated with periodic orbits 

The foregoing reso nances are not the strongest, however. Indeed, we have shown in 
Rieutord et all (|2001l ). that there exist a finite number of frequencies of the spectrum 
which are associated with strictly periodic orbits of the characteristics. These frequencies 
are Up^q = sin(p7r/2g) where p and q are integers; they are such that the angle between 
characteristics and the rotation axis is a rational fraction of tt. Not all the rationals are 
allowed however, because character istics must propagate periodically inside or outside 
the "shadow" of the inner shell f see IRieutord et alll200l[ ). This conditions imposes that 



r] < uj < \/\ — rf for the simplest periodic orbits (it is more restrictive for more com- 
plex orbits). In the case chosen here, namely 77 = 0.35, the allowed periodic orbits are 
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i9 = 


7r/6 


i9 = 7r/4 




Tv/3 




UJ2 


Tl 


Tl 




Tl 


Mode 1 


1.610** 


-2.810^ 


-1.05 10^ 


-210^ 


-2.8 10^ 


Mode 2 


5.8 10^ 


-1.510^ 


-4.32 10^ 


-810^ 


-1.510* 


Mode 3 


3.810" 


-3.810* 


-9.6410^ 


-6101" 


-3.710* 



Table 1 . Eigenvalues of the first modes associated witli periodic orbits of the torus and verifying 
the same symmetry as the forcing (|2.3|l . Each eigenvalue is written A = isini9 + iuj2E^ + tiE, 
as expected for regular modes with stress-free boundary conditions when ui = (the reason for 
this vanishing term is not clear). The ciJ2-term for the 7r/4-modes could not be evaluated because 
of round-off errors. The coefficients n and lj2 have been evaluated numerically using values of 
E around 10"*^. 



associated with the three angles: 7r/6, tt/4 and tt/3. For the associated frequencies the 
trajectories of characteristics are strictly periodic, thus no shear layer is generated (the 
mapping has no focusing power). This situation is illustrated in figure [T] As shown, no 
small scale comes in, and the dissipation diverges as E~-^ when E ^ 0, while the width 
of the resonances diminishes as E. Inspecting the spectral content of the flow, we note 
that the critical latitude contributes to some parts of the flow but at such a low level 

(less than 10"'*) that it does not influence the resonance. 

This feature of the spectrum may be understood as follows. Maas fc Lam ( 19951 ) showed 



that the two-dimensional semi-elliptic basin owns a denumerable set of regular eigen- 
modes. These eigenmodes are associated with periodic orbits of characteristics. They 
also showed that the eigenfrequencies are infinitely degenerate. This is because of the 
ill-posed nature of the eigenvalue problem: For each eigenvalue, the eigenmode is defined 
by an arbitra ry function which i s given on the so-called fundamental intervals of the 



boundary (see lMaas fc Lamlll995l for details). In the case of our fluid domain, the pres- 
ence of an inner core removes almost all the periodic orbits, letting only a finite number 
of them, depending on the size of the core. However, for each of these frequencies which 
are of the form sin(p7r/g) (see above), the associated eigenfunctions are still infinitely 
degenerate since they are specified by an arbitrary function. We illustrate this property 
by computing the first eigenmodes associated with these orbits (see tableland figure [5]). 

The response of the fluid to a time-periodic forcing near these frequencies thus results 
essentially from the superposition of these eigenmodes according to their projection on 
the exciting body force. As the tidal forcing is on the large scale, the responding flow is 
essentially on the large scales. 

2.5.4. Dissipation at resonances 

The foregoing resonances, either associated with attractors (section [2321) '^^ with peri- 
odic orbits (section r2.5.3p . may be described by a simple model of a resonating eigenmode. 
Although much simplified, this model is useful to understand the origin of the scaling 
laws verified by the viscous dissipation. 

Let C be the linear operator governing the forced flow; we write 

iuju = £{u) + f (2.5) 

where u is the frequency of the forcing. Let {uji}„giN be a set of eigenfunctions of C 
verifying the same boundary conditions as u. Hence 



XnUn = £{Un) 
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Figure 7. Left: Rescaled dissipation as a function of the forcing frequency, also rescaled, near 
the periodic orbits with 9ci ~ 7r/3 or uj = \/3/2. Right; the kinetic energy distribution of the 
corresponding forced fiow at i? = 10~^. 
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Figure 8. Kinetic energy distribution of the two first modes associated with the i9 = 7r/4 

periodic orbit given in table [T] 

where A„ = zu;„ + t„ is the associated eigenvalue. We assume that the {un}n^^ form a 
complete orth ogonal basis. As far as modes associated with an attractor are concerned, we 
have shown in Rieutord et ^ ([2002), that they may be described by a Hermite function. 
They thus form a complete basis for the ID functions defined along a line orthogonal to 
the attractor. Thus we may write: 



Here, we assumed that / belongs to the same function space as u. From (|2.5p . we easily 
find that a„ = fn/{iLU — A„) so that 



E 



fn^n 

iuj — Xj. 



(2.6) 
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Now assuming that a single eigenmode dominates the series, we simplify (12. 6p as 

fn^n 
U ~ — 

lUl — An 

The dissipation rate can now be evaluated from (|2.4p . The width of the shear layers is 
assumed to scale like E". The volume that contains the shear layers also scales like E" 
and the gradients of the eigenmode Un scale like E~°'. We get: 

rr r j^^—oi 
D= \c\^dV^- ^— ^ (2.7) 

2 J(V) ~ ^nV + 

Assuming that the complex eigenfrequency of the excited mode expresses as 
A„ — iojn + ME^, we find that 

at the resonance. For reso nances associat e d wit h attractors with vanishing Lyapunov 
exponents, the solutions of iRieutord et al. ( 2002 ) give a = 1/4, P — 1/2 and thus D 



-E"^/^ as observed numerically. For the sharp resonances associated with the periodic 
orbits of section I2.5.3i no small scales comes in so that a = 0. The damping rate of a 
mode with a typical scale independent of viscosity is proportional to E; thus (3 — 1 and 
we get D ^ E^^ as also observed. 

Although not fully rigorous, this short analysis shows that the observed resonances of 
the two-dimensional model behave in a standard way. 

2.6. The critical latitude singularity 

Before ending this section we wish to discuss shortly the role played by the critical 
latitude singularity. We recall that this singularity comes from the " oblique nature" 



of th e boundary conditions to be used with the Poincarc equation (see IRieutord et al 



200l[ ). It leads to a singularity of the solutions but weaker than the one associated with 



the attractors. When the problem takes into account the fluid's viscosity, this singularity 
manifests itself as a broadening of the Ekman layer, which thicke ns from the usual E^^'^ 
scale to the E^^^ scale, on a latitudinal extension that is 0{E^/^) ( Roberts fc StewartsonI 



19631 ). As illustrated in figure H this singularity generates its own network of shear layers 



thus adding some dissipation to the one of the attractor. However, the contribution of 
this singularity is much smaller than the one of the attractor. We also observe that if 
the frequency of the forcing is not such that the attractor has a branch grazing at the 
critical latitude, the amplitude of the flow in this region vanishes with a vanishing Ekman 
number. 

2.7. Summarising the two-dimensional case 

The preceding results show that the dissipation associated with periodically forced shear 
layers may vary very strongly as a function of the frequency of the forcing. At a generic 
frequency, we fi nd that the di ssipation is independent of the viscosity and thus confirm 
the analysis of Qgilvie ( 2005[ ). However, the inertial frequency band contains also in- 



finitely many frequencies (at accumulation points) where attractors are weaker (their 
Lyapunov exponent vanishes and the convergence of characteristics is algebraic instead 
of exponential). At these points the dissipation strongly depends on viscosity, namely 
as E-^/'^. Finally, we also exhibited resonances that are associated with the few allowed 
periodic orbits of characteristics. These resonances are th ose of the regular mo des which 
remain of the dense set of modes of the r/ = torus (see Maas fc Lam 1995f) . As these 



periodic orbits, there is only a finite number of such resonances for a given set-up with 
77^0. 
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We now examine the three-dimensional case so as to determine which of these proper- 
ties remain in this more reahstic case. 



3. Resonances in the spherical shell 

3.1. Numerics 

Turning to the spherical shell problem, we now solve (|1.1|) in spherical geometry. We 
discretize the unknowns and the equations using an expansion of the fields on the spher- 
ical harmonics Y™(9,ip) for the horizontal part and using the Chebyshev polynomials 
on the Gauss-Lobatto colloca tions nodes for the radial part. Details may be found in 
Rieutord &: Valdettarol (1997). We just recall here that the velocity field is expanded as 



+ 00 +/ 



with 



1=0 m=-l 



Op 



TV 



V X m 



where gradients are taken on the unit sphere. Using the same expansion for the body 
force (jl.2p . we find 



12 



^3 



Following iRieutordI ( 19871 ). we derive the equations for the radial functions and 
wf„ (r) from the equations of vorticity and continuity. They read 



A 



i+ir 



dr 



El\tAi{ru'^) - iujAiiru'^) = 



where axisymmetry has been assumed. We have used 
1 



(3.1) 



1 + 1) 

tT 



i^AP - 1 ' ' ^ ' rdr-^' 

The set of equations (13. 1|) is completed by stress-free boundary conditions, which read: 



d ( 
^ t m 

dr \ r 



= 



for the radial functions taken at r = 77 or r = 1. 

In figure [HI we give an example of the spectral content of the forced fiow shown in 
figure [Ml at bottom centre, for E — 210"^°. The spectra show that the truncation 
error of the solutions are O(10~'*), while we estimate d the round-off error (w ith double 
precision arithmetics) to a lower value, in this case (see lValdettaro et a/.ll2007l for a more 
thorough discussion of these error matters). 
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Figure 9. Spectral content of the solution shown in figure [141 middle bottom. Left: Maximum 
Chebyshev spectrum of the velocity field: for a given order of the Chebyshev coefficient we 
plot the maximum absolute value of the coefficient obtained over all the spherical harmonic 
components. Right: The maximum spherical harmonic spectrum of the velocity field: for each 
spherical harmonic order £ we plot the largest absolute value of all the Chebyshev coefficients; 
w refers to while v refers to ru^ (see also .Ricutord fc Valdettara,1997i ). 



3.2. Overall view of the resonance spectrum 

As for the two dimensional case, we first scan the whole inertial band, computing the 
viscous dissipation. The result is plotted in figure [TOl There we note that the response 
curve is very spiky as in the two-dimensional case, revealing very strong variations (six 
orders of magnitude at E=10~^). We plotted the quantity uj^D{u!) to remove the l/cu^- 
divergence at low frequencies, since in this range, the flow is dominated by the velocity 
field u = Ar^ sm6/uje^. The broad shape of the curve does not show a marked symmetry 
with respect to "d = 7r/4 unlike its two-dimensional counter-part (figure [2]). This is ex- 
pected since the curvature terms of the SD-operators contribute to break this symmetry. 
Another striking difference is the absence of resonances at the three frequencies corre- 
sponding to strictly periodic orbits of characteristics, namely •& = 7r/6, 7r/4, 7r/3. On the 
contrary, there is an anti-resonance phenomena. 

Before discussing in more details some specific features of this curve, let us focus on 
figure[TlJ Here, following the idea of ,Ogilvie (2009), we have plotted the dissipation curve 
when a frictional force ~ju replaces the viscous force EAU. As shown by figure [TTl the 
two curves are very similar, demonstrating that the fluid response is mainly governed by 
the underlying Poincare operator (i.e. the operator governing the inviscid problem). The 
meaning of this similarity and the relation with the Poincare operator may be enlighted 
when using this simple frictional force. In this case the flow verifles: 



(iu) + j)u + e 2 X M = — Vp + f 



V • w = 



(3.2) 



Setting A = iw-f 7, this problem, which is completed by the boundary condition il-n — 0, 
may be symbolically written: 
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Figure 10. Scan of the whole inertial band for the viscous dissipation in a spherical shell. 
Vertical dotted lines mark the frequencies allowing strictly periodic orbits of characteristics 
-d = 7r/6, 7r/4, 7r/3. 



'(•:)(;>(T :)(;)•(:) 

or 

{XJ-£p)X 

Here, X = \u,p) and £p symbolizes the Poincare operator. Formally, the solution of 
the velocity field may be expressed as: 

u={XJ-jCp)-^f 
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Figure 11. Scan of the whole inertial band for the dissipation using the viscous force (solid 
line) or the friction force (dotted line). Note the similarity of the curves. 



where (A J — Cp)~^ is a kind of resolvent of the Poincare operator, restricted to a vector 
field, since the operator J is "close" to the identity. The dissipation associated with the 
flow is proportional to the norm of u; indeed, it is 

D = j [ \u\^dV 

J(V) 

Hence, using a restriction of the norm of X to the velocity field, we may write the 
dissipation as 

D^\\{\J-Lp)-^f\\^ 

Now, let us consider the e-pseudospectrum of the Poincare operator. This quantity is 
inde ed very appropriate to dea l with a non-normal operator such as the Poincare one 
(see iTrefethen fc Embred[2005l ) . It is defined as the set of complex numbers A such that 

||(Aj-/:p)-i||>iA 

where we recall that the norm of a bounded linear operator L is the number 

/ \\LX\\\ 
max — - — 

ovcrX y \\X\\ ) 

We therefore see that all the regions of the frequency axis where Z? > belong to the 
e-pseudospectrum of the Poincare operator. Indeed, if D > ^ / then 

l/6<||(AJ-£p)-Vll<ll(AJ-/:p)-i|| 



where we assumed that ||/|| = 1. 

Thus, the dissipation curve obtained with the frictional force gives a partial (one- 
dimensional) view of the e-pseudospectrum of the Poincare operator; computing the 
curve for many 7's would give a view of the two-dimensional subsets of the complex 
plane, which belong to this pseudospectrum. 

Now, the similarity of the two dissipation curves associated with the two damping 
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forces, comes from the nature of the solutions. When these are in the form of a well 
defined attractor, the width of the shear layers is a small scale that singles out so that 
EAu ^ -Ek'^u, where k - E^^^^. 

To conclude this point, we see that the solution of the forced problem show another side 
of the Poincare operator, namely its pseudospectrum, and this quantity is independent 
of both the forcing and the frictional force. The peaks of the dissipation curves offer a 
partial (ID-) view of this quantity. 

3.3. Anti-resonances at periodic orbits of characteristics 

We designate by anti-resonances of the dissipation curve, the frequencies for which the 
dissipation vanishes with the Ekman number. 

We consider the three strictly periodic orbits authorised in the shell with rj = 0.35, 
namely those orbits for which c haracteristics remain either inside or outside the shadow 



of the inner core (see lRieutord et al 2001.) . These are associated with d = 7r/6, 7r/4, 7r/3. 



We plot in figure [T2] the dissipation for these three values and note that it decreases 
following quite closely the law D oc E"^!^ . Such a scaling is obviously reminiscent of 
the critical latitude boundary layer. Actually, as shown by figure fT2l (right), the fluid's 
oscillation is confined along the characteristics emitted by the critical latitude boundary 
layer. 

One may retrieve the scaling of dissipation if we note that the volume of the boundary 
layer surrounding the critical latitude singularity is C'(ii^^/^+^/^) and that the velocity 
field scales as E~^l^ . This scaling appears if we observe that the singular inviscid field 
has a finite normal velocity at the boundary and an infinite tangential velocity (see 
Rieutord et aLllioOlf l. Assuming that u^. is of order unity (as the forcing), boundary 
layers relations imply that the tangential velocity scales as E~^^^ . Using the second 
expression of the dissipation, D = \^ys^Re(u* ■ f)dV, we recover the scaling law that 
we observe numerically. We note that this regularised singularity, propagates along the 
(periodic path of) characteristics, but without any focusing, the resulting shear layer 
widens slowly as we move away from the critical latitude of the inner bounding sphere. 

3.4. Dissipation in the frequency bands with strong attractors 

One of the main results established in two dimensions is that a forced flow oscillating at 
the frequency of an attractor dissipates energy independently of viscosity, provided it is 
low enough. As shown in figure fTSl which gives a zoom on the dissipation curve, we clearly 
see a region of the spectrum (below uj — 0.624) where the dissipation is independent of 
the Ekman number. The curve in figure [T^] (diamonds) confirms this convergence to a 
finite dissipation as viscosity vanishes. 

There is however a surprise in this case. If we look at the Lyapunov exponent in 
the same frequency band we note that something should h appen at the fr e quenc y luq 



where this exponent vanishes. According to the results of iRieutord et all (|200lf ). we 
should expect a resonance, as in two-dimensions. No such resonance occurs, and this 
is independent of the symmetry of the forcing. Inspection of the associated flow (see 
figure [13] right) shows that the attractor is not excited. Rather, the critical latitude 
singularity excites a shear layer, which propagates towards the attractor. This is all the 
more surprising that in this frequency band ther e are many eigenmod es well centered on 



the attractor path (an example may be found in IRieutord et a/.ll200ir ). 

This situation is not unique: it is also the case when we consider the simple attractor 
occupying the frequency band [0.529,0.555]. As shown in figure [Ml no resonance occurs 
there. Actually, if we compare the least-damped eigenmode associated with this attractor 
and the forced flow of the same frequency, which are shown in the same figure, we clearly 
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Figure 12. Left: Dissipation as a function of the Ekman number for four given frequencies. 
'Diamonds': \D{E) — D{E = 10~'')| for a forcing frequency exciting the attractor shown in 
figure 111 and 1 131 at lu — 0.621. 'Pluses', 'Stars' and 'Triangles' show the dissipation at frequencies 
respectively: sin(7r/6), sin(7r/4) and sin(7r/3); the dotted straight lines emphasize the power law 
D(E) oc E^^^. Right: the kinetic energy distribution for a forcing at lj = sin(7r/4) and E=2 10"". 
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Figure 13. Left (top): Zoom of the dissipation in the spherical shell for uj G [0.62,0.64]. Left 
(bottom) : The associated Lyapunov exponents in the same frequency range. Right: The kinetic 
energy of the forced flow at u = 0.621 viewed in a meridional plane. The white line marks the 
path of the attractor. 

see that the attractor is not excited. There too, we see that the shear layer eniitted 
by the critical latitude is strongly excited. It propagates towards the attractor. The 
dissipation curves computed at various Ekman numbers shows that for these values of E, 
no asymptotic regime is reached. Actually, a calculation for a specific frequency shows 
that an asymptotic regime exists but at a very low Ekman number {^10^^^)- We also 
note that the Ekman number below which the viscous dissipation is constant depends 
on the "distance" between the attractor and the characteristics emitted at the critical 
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Figure 14. Top: Zoom of the dissipation curve for various Ekman numbers for uj € [0.5,0.6]. 
The two vertical solid lines delineate the frequency range of the attractor shown by the eigen- 
mode below. The dotted vertical line shows the frequency of the least-damped eigenmode. Be- 
low left: The kinetic energy distribution in a meridional plane of the least-damped eigenmode 
{uj — 0.554838) associated with the attractor governing this frequency range. Below centre: The 
amplitude of the forced flow at the frequency of the eigenmode. Below right: The frequency of 
the forcing is now lo — 0.542 so that the attractor is closer to the critical latitude. Note that a 
second attractor is now visible and excited by the southern branch of the shear layer emitted at 
the critical latitude. 

latitude. We observe that at = 0.542 the asymptotic regime is almost reached; in this 
case, the attractor is much stronger and closer to the critical latitude (as confirmed by 
figure [ri]). 

Back to the ui = 0.621 case, we note that all dissipation curves at some frequencies 
above lvq = 0.622759 (which is the upper limit of the attractor living at w < loq), seem 
to have converged to a limit. We think that the asymptotic limit has not actually been 
reached there, and that undulations of the dissipation curve will appear, as they do at 
higher frequencies, in response to the rapid variations of the Lyapunov exponent. 

To conclude this point, the foregoing examples show that the asymptotic regime, where 
dissipation is independent of viscosity, is reached when the Ekman number is low enough 
so that the shear layer emitted at the critical latitude on the inner boundary is "feeding" 
an attractor. In this case the widening and softening of the shear layer, due to viscous 
diffusion, is stopped. The shear layers are then in a regime similar to what has been 
described bv iOgilvia (2005.) in two dimensions. 
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Figure 15. Zoom of the dissipation in the spherical shell in a given frequency range; the four 
vertical straight lines show the range of existence of two short-period attractors. 



3.5. Resonances 

The two-dimensional case gave a nice illustration of resonances which correspond to 
frequencies where the Lyapunov exponent vanishes. The foregoing discussion has shown 
that the three-dimensional situation is not so simple. Actually, we could not identify a 
single resonance that matches the frequency corresponding to the vanishing Lyapunov 
exponent of a well-determined attractor! 

Nonetheless, dissipation curves show many peaks indicating some privileged frequen- 
cies. The formation of these peaks as viscosity decreases is not quite standard. Indeed, 
if we compare the dissipation curves in the interval [0.5,0.6] (fig. I14p and the one in the 
interval [0.65, 0.67], (fig. [15]) there is a common feature. This is the fact that "resonant 
peaks" are more and more numerous as the viscosity decreases (as expected), however 
they do not seem to appear because the resonances are more and more vigorous (it is 
only marginally the case), but because the neighbouring frequencies resonate less and 
less. Hence, it seems to be that many peaks are not true resonances, but small intervals 
of frequencies where the dissipation is independent of the viscosity provided it is low 
enough. The series of peaks for uj £ [0.56,0.58] is quite illustrative of this phenomenon. 

However, there are also some peaks which do behave like true resonances, i.e. their 
amplitude steadily increases as the Ekman number decreases. This is for instance the 
case of the two major peaks of the spectrum, lying at uji — 0.6598 and uj2 — 0.6629. 
These resonances corresp ond to the least-d amped axisymmetric inertial modes that we 
actually studied in .Rieutord fc Valdettarj (1997). We therefore extended this previous 
study to lower Ekman numbers. In figure I16[ we focus our attention to the first mode, 
which seems to follow some asymptotic regime (as shown, the damping rate of the second 
mode is clearly not a power law). Figure [TBI shows that the dissipation, there, grows like 
£;-i/3 when 10"^, but like E-^/'^ if £;< lO^^. On the other hand, the width of the 
resonance narrows as E'^^'^ when E<10~^; we note that the associated eigenmode has a 
damping rate decreasing as Applying the simplified model of section [2.5.41 to this 

case, we recover the scaling law of the viscous dissipation at £'>,10^^ if we assume that 
the critical latitude singularity emits a shear layer of width E^^^, i.e. corresponding to the 
latitudinal extension of the perturbed Ekman layer on the inner sphere. Hence, dissipation 
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Figure 16. Top left: The viscous dissipation at the resonance uj = 0.6598 as a function of the 
Ekman number; the straight lines shows the power laws E~^/'^ and E"^''^. Top Right: The solid 
hne and pluses show the width &uj of the resonance compared to the power law E^'^'^ (dotted 
line). Below: We show the damping rate r of the two least-damped eigenmodes at u\ = 0.6598 
(pluses) and tJ2 ~ 0.6629 (soUd line); the power law E"'^'' is shown by the dash-dotted straight 
line. 



should scale as like £:4/5-2/3_ ^gj^g ^ ^ q 57 

we find D ^ E A similar matching 
of the exponents can be obtained if the forced flow is confined to a neighbourhood of 
the critical latitude. The volume responding to the forcing is 0{E^/^) but the velocity 
gradients are 0{E~'^/^), thus leading to the same power law for the dissipation. The 
simple model of section \2 . 5 . 41 cannot be more precise, but it clearly u nderlines the cru cial 
role played by the critical latitude. We refer the reader to the work o f lKerswelll (|l995l) for 
a detailed analysis of the boundary layers and shear layers in the vicinity of the critical 
latitude in the case of the spin-over mode (i.e the non-axisymmetric, m=l, aX lj ~ 0.5, 
inertial mode). 

Now, the power law change at E' ~ 10~^ may be understood with figure [171 In this 
figure we show an enlarged view of the resonance together with the Lyapunov exponents. 
It is clear that, as the Ekman number decreases, the resonance shifts to the frequency 
interval where this exponent is very small. We note that E = 10"^ corresponds to the 
transition where the resonance leaves the frequency range occupied by a short-period at- 
tractor and enters the frequency interval where only long-period attractors exist, leading 
to a stronger resonance. The new i<^~^/^-regime is likely not asymptotic as well: Inspection 
of the Lyapunov curve shows that the size of the frequency intervals, where a definite Lya- 
punov exponent exists, is at least less than 10~*. This means that an asymptotic regime 
may be reached at Ekman numbers less than 10~^^ (assuming that the resonance keeps 
narrowing as E'^/'^). Such numbers, even if theoretically reachable in stars or planets, are 
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Figure 17. Top: The viscous dissipation at the resonance near uj = 0.6596 as a function of the 
frequency for various Ekman numbers. Below: The corresponding Lyapunov exponents in the 
same frequency interval. Note the progressive move of the resonance peak towards the frequencies 
with small Lyapunov exponents. 



likely unrealistic because any small-scale turbulence would increase them by many orders 
of magnitude. Thus, the intermediate regimes revealed by the E~'^/^ or laws for 

the dissipation are likely to be more relevant to astrophysical or geophysical applications, 
but more investigations are needed to determine their origine. 

3.6. Influence of the size of the core 

A last question that is often raised is whether the size of the core, when it is small, 
influences notably the response curve that woul d be obtained i f negl ecting the core. 
From the foregoing results and the discussion of Rieutord et al. (2QQ0), the answer to 



this question obviously depends on the Ekman number. It is indeed expected that for 
large viscosities, a small inner core is hardly seen by the fluid motion. To illustrate 
further this point we plot in figure 1181 the dissipation curves for various sizes of the 
inner core together with the one of the full sphere. As shown, at a rather large Ekman 
number E — 10^^ (unrealistic for planetary or stellar applications), the dissipation curve 
of the spherical shell notably differs from the full sphere case only when r] > 0.2. On 
the other hand, with a more realistic value E = 10~*, we note that a very small core 
rj = 0.05 already amplifies the dissipation by an order of magnitude. One also observes 
that the frequencies of some resonances are quite similar. This likely comes from the 
polynomial nature o f the full spher e solutions. Using a truncated solution in series of 
spherical harmonics, iRieutordI (|l99lh showed that the frequency of one of the large-scale 
modes of the full sphere was shifted by an amount 0{r]^). 

In fact, in the inviscid case, the presence of a small core, although fully perturbating 
the eigenvalue spectrum of the Poincare operator, likely has a less drastic influence on 
the pseudo-spectrum, which we partially view through the dissipation curve. 

To be more quantitative and determine the threshold rjc, for a given E, beyond which 
the core cannot be ignored, we computed the dissipations at various E, keeping uj fixed, 
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Figure 18. The resonance curve for various size of the inner core and two Ekman numbers. 
Left: Note that the pluses correspond to the case r) = 0.1 and are almost sitting on the curve 
= 0; the difference is always less than two percents. Right: The same scan at a much lower, 
but more realistic Ekman number. A small core induces important differences. 
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Figure 19. Variations of the dissipation at a given Ekman number for varying sizes of the 
inner core. The frequency of the forcing has been arbitrarily fixed to lo = 0.55. 



and varying the size of the inner core. The curves in figure [T9l illustrate these variations 
of the dissipation. They have been computed for lo = 0.55 but the behaviour seems to 
be generic: We picked up two other frequencies and found the same trends. However, 
for u) = yjijl, which is the broad resonance showing up i n figure | 18 l (right ) and which 
corresponds to a large-scale mode of the full sphere (e.g. iGreenspanI 
find this behaviour. 

From the curves in figure [1^1 we could determine that for sufficiently small jy's 



1969I) . we do not 
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D{j^)^D{r] = 0) + a{T^)E + bTt' (3.3) 

where D{ni) is the viscous dissipation, 6 is a constant and a(r/) is a rapidly varying 
piecewise bounded hnear function. The -E-dependence is expected from a flow which 
does not (or little) depend on the viscosity as it is the case for the full sphere m odes. 
Thejj^-dependence can be, tentatively, explained using the solution derived in RicutordI 
( 19911 ) (see the discussion in appendix C). From the law (|3.3p we therefore conclude that 
the inner core influences the dissipation when hif' > OmE, a^n being the upper bound of 
0(77), showing that the critical radius of the inner core is 

Vc = qE^'^ (3.4) 
where the coefficient q weakly depends on the frequency of the forcing. We find q ^ 
0.4. This formula shows that only extremely small cores can be neglected in realistic 
conditions where E<10^^. 

3.7. Discussion 

The foregoing results show that the three-dimensional flows are significantly different 
from their two-dimensional counter part. The numerical solutions show a kind of exchange 
of the roles between the attractors and the critical latitude singularities as one shifts from 
the two-dimensional to the three-dimensional problem. In two dimensions the attractors 
are clearly dominating the dynamics and control the periodically forced flows. The shear 
layer emitted at the critical latitude is present but not essential. In three dimensions we 
observe the opposite: the shear layer emitted at the critical latitude on the inner sphere 
plays a crucial role in the fluid's response to the periodic forcing, while attractors appear 
when they are "fed" by this shear layer. 

We have no definite explanation to this observation. We conjecture that it is a con- 
sequence of two modifications in the solutions when one changes the dimension of the 
problem. The first may come from the different Ricmann functions of the two- and three- 
dimensional problems. We recall that the general solution of an hyperbolic problem may 
be expressed as 



n(^) = i(n(P) + n(Q)) 



+ i / R{S, M) (^du+ - ^du\ + n(M) (^du+ - ^du\ (3.5) 

where R{S; M) is the Riemann function associated with the Poincare operator; inte- 
gration is done on the data line PQ. In 2D this function is a combina tion of a Dirac 



distribution and the identity, while in 3D it is a Legcndrc function (see iRieutord et 



[2001.') . Hence, in 2D, a singularity on the data line will remain on the web of characteris- 
tics issued from this point, while in 3D it may generate singular points outside this web. 
Thus, the influence of the critical latitude singularity in the three-dimensional problem 
is likely stronger than in the two-dimensional one. 

The second modification brought about by the third dimension is the boundary of 
the domain. While in 2D the true fluid domain is a disk with a core, which is not a 
simply connected fluid domain, in 3D the meridional section of the fluid layer is a simply 
connected domain. The difference comes from the presence of the rotation axis which is a 
line of the meridional plane where the characteristics must reflect, which is not the case 
in 2D. We conjecture that this difference makes the eventual existence of fundamental 
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intervals (introduced by iMaas fc Lam 1995 ). more difficult thus making the excitation 
of shear layers along the attr actors also more difficult. 



4. Conclusions 

To conclude this work, we may summarise our results in the following way. First, it 
turns out that the two-dimensional model, which is equivalent to the section of a slender 
(cored) torus, is well understood: 

• The dissipation curve is very spiky as one scans the frequency axis and each spike 
corresponds to a resonance. This resonance is associated with a frequency for which the 
Lyapunov exponent is zero. We have shown that in most cases the flow is associated with 
an attractor whose shear layer width scales as E^/'^. In a few cases, when the frequency of 
the forcing is of the form sin(p7r/2q), where p and q are integers, the resonance is stronger, 
but narrower, and corresponds to one of the few eigenmodes which remain when an inner 
core is inserted at the centre of a circular domain. 

• In between resonances, which are not dense (but have accumulati on points), th e 
dissipation reaches a constant value at vanishing viscosity, as predicted bv lOgilvi^ ( 2005 ). 



The two-dimensional model thus offers a rather neat picture well constrained by the 
analytical results. It shows that, contrary to previous expectations, dissipation is very 
sensitive to viscosity, essentially because of resonances but also because, outside reso- 
nances, shear layers associated with different attractors reach their asymptotic regime at 
very different Ekman numbers. 

In three-dimensions, results show some similarities with the two-dimensional model 
but also major differences. Among the similarities, we note that the dissipation curve is 
very spiky too. There also exist intervals of frequencies where the dissipation does not 
depend on viscosity when this quantity is low enough. There are also resonances where 
the dissipation seems to increase without bounds when the viscosity vanishes. 

However, there are major differences. First, we observed that the critical latitude sin- 
gularity on the inner sphere plays a major role in the response of the fluid to the forcing. 
It is systematically emitting a shear layer. Using the scalings of the Ekman layer at this 
place, we could explain the vanishing dissipation with vanishing viscosity when the fre- 
quency of the forcing is associated with a periodic orbit of characteristics. The latter case 
also stresses another difference: whereas in the two-dimensional case these periodic orbits 
were associated with strong and narrow resonances, in the three-dimensional case they 
are associated with anti-resonances : the dissipation vanishes with the viscosity (following 
a i;2/5-law). 

We observed that in the two-dimensional case, least-damped modes are associated with 
frequencies where the Lyapunov exponent is zero. In three dimensions, such an associa- 
tion also exists but some of the least-damped modes are also found in association with 
intervals of frequencies where the mapping is weakly contracting (because of long-period 
attractors). Our numerical results on the forced problem show that only this latter kind 
of least-damped modes lead to strong resonances. The question of whether these reso- 
nances can reach an asymptotic regime remains open; the fractal nature of the Lyapunov 
exponents in the frequency region where they exist may prevent any asymptotic limit. 
Clearly, for one of them, even a,t E = 10~^^, the asymptotic regime is not reached. 
The close examination of this resonance has nevertheless shown that some intermediate 
regimes may exist when attractors of very long period cover the same intervals of fre- 
quencies as the resonances. Such regimes are astrophysically relevant and deserve more 
investigations. 
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Unfortunately and unlike in the two-dimensional case, we have little analytical guide 
in the three-dimensional case. Tentatively, we interpreted the striking role of the charac- 
teristics emitted at the critical latitude as a consequence of the nature of the Riemann 
function associated with the Poincare operator and the probable different fundamental 
intervals associated with the mapping of characteristics. 

Finally, the striking general conclusion is that a periodically forced flow in a spherical 
shell is quite different from what can be expected from a simple response of resonant 
eigenmodes. Previous work has shown that eigenmodes were often featured by attractors 
of characteristics, but those modes turned out to be of little interest for the forced prob- 
lem. Interestingly enough, the dissipation curve reveals some parts of the pseudospectrum 
of the Poincare operator, thus offering another way to investigate the properties of this 
operator. 

The solutions that we computed also showed that the critical latitude singularity is a 
determinant feature of the periodically forced flows. More work is now needed to fully 
understand the interplay of the singularities ge nerated by both the cr i tical l atitude and 
the attractors of characteristics. We note that Goodman fc Lackner (|2009 V while ex- 



amining the tidal response of a Jovian planet, were also led to the conclusion that the 
critical latitude at the interface between a solid core and a fluid envelope plays a crucial 
role in the dissipation of kinetic energy. In this same planetary context, we showed that 
planets with a very small core could be assimilated to a full sphere only when the relative 
radius of the inner core, rj, is less than 0.4_E^/^. Beyond this, usually very small, value 
the dissipation grows rapidly with the core size as if . 

To conclude on the astrophysical problem, which motivated this study, our results 
stress the importanc e of viscosity. They unfortunately remove the possibility, envisaged by 
Qgilvie fc Lin liooi), that the synchronisation time scale of binary stars be independent 



of the viscosity, which is usually not a well known quantity. 
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Appendix A. Viscous dissipation in the slender torus 

A.l. Equivalence of coordinate systems 
Using non-dimensional variables, viscous dissipation is 

D = ^ I \c?dV 



2 



{V) 



where |c| symbolises the norm of the shear tensor. In usual cylindrical coordinates, this 
norm reads 

\C? = CijC*j = \Css\^ + \c^^\^ + \c^^f + 2(|c,^|2 -h \c,^f + |c^^|2) 

with 

. -2^ c -2/^i^ + ^^ c -2^ 
as \s dip s J oz 

_ldv^ dv^^v^ _ dv^ dvs_ _ dv^ 1 dv^ 

s dip OS s OS az oz s dip 

Considering axisymmetric solutions aA, 1, the shear tensor components simplify into 

-9^ -n _ c^dvz 

_ dv^ _ dvz dvs _ dv^ 

ds ds dz dz 

Using the meridional stream function ijj, we may further write 

^zz — ^ a a — ^ss 

dsdz 

_ d^tjj d'^tp 
ds'^ dz^ 

du du 

= -g^, C^Z = 

So that the volumic dissipation is : 



Now a direct calculation shows that: 
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and 



where we recognise the expressions of Cpp and Cp^ of (jcj). Thus, 

c^^ = -Css = sin2(/>Cp0 + cos20Cpp 
= cos 20Cp0 - sin 2(t>Cpp 

In the same way, we also find that 

= cos 4>Cpi^ — sin 4>c^c, and c^^ = sin (jyCpQ + cos (/>C0^ 

which leads to 

|cp = 2(|cppp + |cpcP + |c^Cp + |cp^n 

as expected. 

A. 2. Dissipation as function of the Fourier components 
We first note that 



/ |c,cP + |c,cPd0 = 27r5^|K:P + ^|K 

r \cpp\'de = 2nJ2 
/ |cpepd^ = 27rV 



2n / , , ■0« 

— h/'n 



p p 



Since the total dissipation is 



D = E (Icppp + IcpcP + \cec\' + \cpe\^)dV, 

J{V) 



it finally expresses as 
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— hAn 

P \ P 



P P 

If the solutions are equatorially symmetric — Vn and V'-n = ""^n: so that 



pdp 



-'^ n>0 P 



Wn 



— hZ-n 

P \ P 
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P P 



Appendix B. Influence of boundary conditions 

iFotheringham Sz HollerbachI ( 1998h studied the infiuence of boundary conditions (no- 
slip or stress- free) on the inertial modes of a spherical shell. They found that this infiuence 
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was very small. In the forced flow that we study here, this is also the case. The contribu- 
tion to dissipation of Ekman layers that appear if no-slip boundary conditions are used 
may be estimated as follows: In the case of the asymp totic regime w here the fluid flow 
is following an attractor, the velocity scales as E~^/^ [ Qgilviel 200 5^ and the boundary 
layer velocity gradients du/dx are 0{E-^'^ /E^/'^) = 0{E-^/^). However, the volume 
occupied by the boundary layers is 0{E^/^)^ thus, the resulting dissipation scales like 
is very small compared to the contribution of the internal shear layers, which 
is unity. In the case of the resonances of the 2D-slendcr torus, which are associated with 
shear layers of width scaling in E^/^^ the contribution of Ekman layers is even weaker, 
namely 0{E^/^). It is only in the case of the regular modes of the torus, that a strong 
effect is noticed, as expected. However these modes do not exist in the more realistic set- 
up of a spherical shell: we obtained singular solution (the "anti-resonances" see § I3.3p . 
In this case we find that the dissipation scales like E"^/^ , which seems to be also the case 
when a no-slip boundary condition is used on the inner core. 



Appendix C. The asymptotic law for cores of vanishing sizes 

We may recover the depend ence in rf' of the viscous dissipation when 77 <C 1 by con- 
sidering the solutions found by iRieutordI (|l99lh . Indeed, in this paper it was shown that 
an oscillatory flow in a rotating background with spherical symmetry could be described 
by two kinds of solutions when these are decomposed onto the spherical harmonic basis. 
These solutions (in fact the radial functions of each spherical harmonic component) are 
either of Bessel or polynomial type. The Bessel type solutions describe the boundary 
layer regions, while the polynomials describe the flow in the remaining volume^ 

Con sidering the inviscid case, we shall focus on the polynomial solutions. iRieutordI 
(|l99l[ ) has shown that these solutions also split into two categories: one is regular at the 
origin, the other is regular at infinity. The first one is e asily related to t he inertial modes 
of the full sph ere which have been found bv iBrvan ( 1889^ (see also Greenspan 19691 : 
Rieutordl[T997h . They are exact solutions of the inviscid problem at eigenfrequencies. 



Below, we shall refer to these solutions as class-one solutions. 

The second class of solutions are those with radial functions with polynomials in 1/r 
that are regular at infinity. However, unlike those regular at the origin, these solutions 
are exact only when the momentum equation is expanded on a finite number of spherical 
harmonics, namely when the series is artificially truncated at a given maximum order, just 
like in a numerical calculation. This is expected since the solutions regular at infinity ar e 
singular at the critical latitude on an inner bounding sphere (see iRieutord et a/.l 120011 ) . 
Below, we refer to these solutions as class-two solutions. 

Now if we think to the solutions in their spherical harmonic decomposition, we note 
that the radial velocity of the class-one solution is dominated near the centre by its u^l2 
component (we restrict the discussion to the axisymmetric and equatorially symmetric 
flows). This is because ~ r^~^ as r 0. When a small inner core is introduced in 
the inviscid problem, we need to add solutions that are regular at infinity so that the 
inner boundary condition Ur{r]) = is met. This is the usual procedure in regular elliptic 
problems like the Poisson problem for instance. However, in the Poincare problem, this 
cannot be done because the class-two solutions are singular. However, if we restrict our 
problem to a finite number of spherical harmonics (which is a way of regularization!), we 
can use this ki nd of solut i ons. 

As shown in IRieutordI (1991), the form of class-2 solutions for w^(r) is in r~'*. If we 
write 
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u = Uq + a2uio + oii^to + • • • 

where uq is the class-1 solution and are the class-2 solution of order J (their lowest 
order is (r) = r~'^). Imposing that Ur{ri) = for r/ <C 1, we easily find that a2 is 0{ri^) 
and that more generally, a2p is 0(77^^+'^). 
Hence, to leading order, we may write 

u — uo + ifun (C 1) 

Although the truncated solutions are not exact, we know that viscosity will, in the end, 
also truncate the spherical harmonics series. Computing the viscous dissipation from 
(|C ip shows that 

D{'n)^D{'n = Q)+rfDi 

where the Z)i-term might depend on the Ekman number. Although quite close to the 
searched expression (|3.3|) . this equation misses the a{rj)E term and the fact that Di 
seems to be independent of E. 

We note that the critical latitude boundary layer on the inner sphere could induce 
a velocity field which is 0{r]/E^^^). This scaling also suggests that a transition occurs 
when 77 ^ E^/^. This latter remark shows that the complete explanation is likely involved 
and a more detailed analysis is necessary to fully understand the intricacy of the limits 
7y ^ and £; ^ 0. 



